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ABSTRACT 

As a first step towards understanding the angular momentum evolution history of black holes in 
merging black-hole/neutron-star binaries, we perform population synthesis calculations to track the 
distribution of accretion histories of compact objects in such binaries. We find that there are three 
distinct processes which can possibly contribute to the black-hole spin magnitude: a birth spin for 
the black hole, imparted at either (i) the collapse of a massive progenitor star to a black hole or (ii) 
the accretion-induced collapse of a neutron star to a black hole; and (iii) an accretion spin-up when 
the already formed black hole [via (i) or (ii)] goes through an accretion episode (through an accretion 
disk or a common-envelope phase). Our results show that, with regard to accretion-induced spinup 
in merging BH-NS binaries [method (iii) above], only accretion episodes associated with common- 
envelope phases and hypercritical accretion rates occur in the formation history of merging black 
hole/neutron star binaries. Lacking unambiguous experimental information about BH birth spins 
[i.e., regarding the results of processes (i) and (ii)], we choose two fiducial values for the BH birth 
angular momentum parameter a = consistent with observations of (i) NS birth spins (a ~ 0) 

and (ii) X-ray binaries (a = 0.5). Using these two fiducial values and a conservative upper bound 
on the specific angular momentum of accreted matter, we discuss the expected range of black hole 
spins in the binaries of interest. We conclude with comments on the significance of these results for 
ground-based gravitational-wave searches of inspiral signals from black hole binaries. 

Subject headings: binaries:close — starsievolution — stars:neutron - black hole physics - gravitational 
waves - relativity - X-rays:binaries 


1. INTRODUCTION 

Gravitational-wave astronomy is a rapidly-developing 
field, with many ground-based interferometric detectors 
(LIGO, VIRGO, GEO, TAMA) in place and taking data. 
These detectors are sensitive to (among other things) the 
waves emitted from the late stages of inspiral and merger 
of compact binaries with neutron stars (NS) and black 
holes (BH). Under certain conditions (as discussed in 

A [jOsto^os_et al.lll994 iGrandclement et 'mi200.‘li 12004 

Pan et musir and references therein) the BH spin in 

high-mass ratio binaries such as BH-NS binaries can 
significantly modulate the emitted inspiral waveform. 
While the more complex resulting waveforms can pose 
challenges for detection strategies, they also enable us to 
empirically determine the spin of rapidly-rotating black 
holes, if they exist. We undertake the present study of 
the expected BH spin distribution for two reasons: (i) 
to examine whether special search methods for spinning 
compact objects for optimal detection efficiency are nec¬ 
essary from an astrophysical point of view; (ii) to provide 
the first steps for a theoretical understanding of the ori¬ 
gin and magnitude of BH spins in BH-NS binaries that 
will be useful when such BH spins are empirically con¬ 
strained in the future. 

In Sec. El we briefly describe the population synthesis 
calculations we use to explore a broad range of possi¬ 
ble scenarios that lead to the formation of BH-NS bina¬ 
ries that can merge within a Hubble time (hereafter re- 
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ferred to as “merging”). We summarize and explain the 
common formation channels for these binaries and de¬ 
scribe the mass accretion history associated with them. 
In Sec. El we review the observational evidence for BH 
birth spin drawn from observations of isolated NS and of 
X-ray binaries. Based on this evidence, we employ two 
fiducial choices for BH birth spin parameter a = J/IVP, 
i.e., 0 = 0 (nonspinning) or a = 0.5 (moderate spin). 
Given an understanding of the mass accretion history 
and a choice for the initial BH spin , we could in prin¬ 
ciple determine the final black hole spin if we knew the 
specific angular momentum of the material accreted by 
BH in these systems. This factor remains substantially 
uncertain. Thus, in Sec. El we use a conservative model 
for the specific angular momentum, combined with the 
distribution of accretion histories and birth spins, to limit 
the rate at which high-spin black holes occur in merging 
BH-NS binaries. 

2. ACCRETION HISTORY OF BLACK HOLES IN BH-NS 
BINARIES 

To generate and evolve stellar populations until dou¬ 
ble compact-object formation occurs, we use the Star- 
Track code first developed by Belczynski, Kalogera, 
and Bulik (2002) [hereafter BKB] and recently signif- 
icantly updated and t ested as described in detail in 
iBelczvnski et all ll2005l . As with all binary popula¬ 
tion synthesis codes, StarTrack encapsulates the residual 
physical uncertainty in evolutionary processes of stars 
in a number of parameters, several of which can signif¬ 
icantly affect the statistics of BH-NS mergers. In order 
to have a sample of merging BH-NS drawn from a large 
k&eiJiijY®QfeLa$.tiBSii)hysically a priori plausible models, we 
generated an archive of accumulated results from a broad 
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swath of plausible models for population synthesis. As 
described in O’Shaughnessy, Kalogera, and Belczynski 
(2004) [hereafter denoted 0KB], the archive was gener¬ 
ated using flat probability distributions for seven of the 
most significant parameters: three parameters describ¬ 
ing the supernova kick magnitude distribution (assumed 
to consist of two maxwellian distributions with disper¬ 
sions (Ti S [0,200]km/s and U 2 G [200,1000]km/s and 
relative weight of small to large kicks s G [0,1]); the ef¬ 
fective common-envelope efficiency aX G [0,1]; the stel¬ 
lar wind strength w G [0,1]; the power-law index con¬ 
trolling the distribution of the mass of the companion 
relative to the mass of the primary r G [—3,0]; and the 
fraction of mass accreted (as opposed to lost from the 
binary) during non-conservative mass transfer episodes 
f„ G [0,1]. Unfortunately, as discussed in more detail in 
lO’Shaughnessv et"^ ll2005bll . we cannot succinctly de¬ 
scribe the effect of these seven parameters on the forma¬ 
tion channels and rates of BH-NS systems; physical pre¬ 
dictions depend in a highly correlated way on all seven 
parameters. As in 0KB, we fixed the remaining popula¬ 
tion synthesis model parameters to physically reasonable 
values, such as setting the maximum neutron star mass to 
M]ys,max = 2 Mq and the metallicity to solar metallicity 
Z = 0.02 (see parameters of model A in BKB for more 
details). From this archive we identified each merging 
BH-NS binary along with the details of its evolutionary 
history. Our results are summarized in Tabled 

Compact binaries tight enough to merge through the 
emission of gravitational waves within the Hubble time 
have progenitors which usually interacted strongly (e.g. 
through mass transfer) earlier in their evolution. For 
example, a common mode of interaction is conventional 
mass transfer through Roche-lobe overflow and disk ac¬ 
cretion. However, for merging BH-NS binaries, we And 
in our simulations that almost all of them have expe¬ 
rienced a more dramatic phase of dynamically unstable 
mass transfer and common envelope (CE) evolution (see 
Table 0. During such phases it is expected that hyper¬ 
critical accretion becomes possible. In t his form of ac ¬ 
cretion (HCE), disc ussed in more deta i ] inlBrownl (|1 flflfJl . 
iBrown et al.1 i200(11 . and iFrver et al.1 ll 1 99(111 . a compact 
object spirals in through the envelope of its companion, 
rapidly accreting matter at highly super-Eddington (for 
photons), neutrino-cooled rates. As a result of this pro¬ 
cess, most of the companion’s envelope is ejected, bring¬ 
ing the post-CE binary very close together. While the de¬ 
tailed accretion onto the compact object is complex and 
ill-understood, since we know the whole envelope will be 
lost during this phase, a straightforward application of 
the conservation of energy and mass during the quasicir¬ 
cular spiral-in process can determine the final orbit, the 
final binary masses, and the mass of the eject ed material; 
see Appendix A of iBeIczvnski et ^ l^2^^2^ for details. 
Nevertheless we stress that our treatment of HCE is es¬ 
sentially a simple, semi-analytical model, which is by no 
means guaranteed to accurately describe the quantitative 
effect of the HCE phase. 

More specifically, as described in Tabled the vast ma¬ 
jority (all but 582, out of 5435) of BH-NS binaries seen 
in our simulations form through an evolutionary channel 
that involves an HCE phase. In these channels, after an 
optional mass transfer phase, one star explodes, the com¬ 
pact remnant spirals through and strips the envelope of 


TABLE 1 

Statistics for evolutionary channels 


Channel 

n 

AIC 

(MT-l-)SNa 

HCE(b^a) 

(MT+)SNb 

4310 

3656 

(MT-l-)SNb 

HCE(a^b) 

(MT+)SNa 

543 

543 

(MT+)SNx 


MT+SNy 

319 

0 

(MT+)SNx 


SNy 

263 

0 

total: 

5435 


Note. — This table summarizes all evolutionary channels fol¬ 
lowed by systems which formed merging BH-NS binaries in our 
simulations. The first three columns describe the evolutionary 
channel (where SNx indicates a supernovae of either a the primary 
or b the secondary, HCE indicates a hypercritical common enve¬ 
lope phase, and MT indicates a stable Roche-lobe-overfiow mass 
transfer; parentheses indicate an optional feature of the channel); 
the fourth column provides the number of merging BH-NS binaries 
which passed through this channel; and the fifth lists the number 
which undergo accretion-induced-collapse (i.e., the BH forms dur¬ 
ing the common-envelope phase from an accreting neutron star). 
Most systems formed via a hypercritical common envelope (HCE) 
phase onto a compact object; this compact object is occasionally 
already BH, but usually is a NS which undergoes accretion-induced 
collapse (AIC) to a BH. A smaller number of systems (582 of 5435) 
do not undergo HCE at all and interact either weakly (through 
winds) or not at all between the two supernovae. 


its companion, and then the companion explodes. Black 
holes form either immediately after the first core-collapse 
event or through the accretion-induced collapse of a NS 
that forms in the first supernova event. Regardless of 
the nature of the first compact object (BH or NS), it 
most typically experiences HCE evolution, and accretes 
a large amount of matter (potentially equal to its birth 
mass) during the brief HCE phase. EigureQ shows the 
fraction of the initial compact object’s mass M-mit that 
is accreted during the HCE phase. Based on the sample 
we have, we estimate an upper bound on the fractional 
mass accreted as a function of the initial mass to be 

(AM/M) max « 0.25 + 2MQ/Minit . (1) 

This bound is valid for all channels.^ Eor systems which 
form their BH immediately after the first supernova (i.e., 
no accretion-induced collapse, so Ma,init > 2Mq) we And 
an empirical lower bound of 

(AM/M) max « 0.05 + O.GMe/Minit . (2) 

This latter bound may be understood as a requirement 
that the BH-NS system merge within the simulation 
time. Given that the second supernova will rarely fur¬ 
ther tighten the orbit, if the binary is to merge within 
10 Gyr, then the orbit must be smaller than roughly 
3.ARq[{Mi/Mq){M2/Mq){Mi + M 2 )/Mq]^^'^ after the 
common envelope phase. Using the common-envelope 
evolution equat i ons p resented in Appendix A of 
IBeIczvnski et alJ (I2002D (with any choice for acE x A), 
we And this requirement approximately translates into 
the lower bound presented above. 

Merging BH-NS can form through other channels, 
as shown in Table Q] However, these channels in¬ 
volve substantially less mass transfer onto the BH: for 
the HCE-related channels, the smallest mass increase 

® We expect that systems which transfer more mass than this 
limit are brought so close that they merge during their common- 
envelope phase. 
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Fig. 1.— For systems that undergo HCE evolution (see Table0, 
this plot shows (Mfinal - MnitaO/Mnitial versus Minitial, where 
Minitial Mfinal t^e mass of the BH progenitor just after its 
supernovae and at the end of its evolution, respectively. In other 
words, they show the relative increase in mass of the compact ob¬ 
ject AMhce/-^ due to hypercritical common envelope accretion. 
The upper solid line shows the empirical upper bound given in 
Eq. Q; the two vertical lines shows our assumed lower and upper 
limits on NS mass, M = I.SMq, 2Mq. 


was 0.22Mq, to be contrasted with the largest disk- 
mediated mass transfer, O.OOTMq (i.e., from the chan¬ 
nel SNx-l-MT-l-SNy). As described in Section 0] this in¬ 
credibly small disk-mediated mass transfer cannot signif¬ 
icantly modify the birth spin of the BH. Therefore, we 
neglect these channels in subsequent discussion. 

Caveat: Branching ratios: Unfortunately, our branch¬ 
ing ratios for various evolutionary channels are not sig¬ 
nificant. When exploring different population synthesis 
parameter combinations, we normalized so each popu¬ 
lation synthesis run contributes a roughly fixed num¬ 
ber of merging NS-NS binaries to the archive. As a re¬ 
sult, models associated with low NS-NS merger rates but 
high BH-NS merger rates should contribute dispropor¬ 
tionately more merging BH-NS to our archive, and vice- 
versa. Moreover, even if each model contributed equally 
to the total number of BH-NS systems, our net num¬ 
bers would then reflect the mean branching ratios over 
all population synthesis models considered, rather than 
the branching ratio o f. e.g., the most common mode l. Fi¬ 
nally, as discussed in lO’Shaughnessv et al ] ll2005bll . not 
all population synthesis models we considered are equally 
compatible with observations. Nonetheless given the cur¬ 
rent statistics, we fully expect that most BH-NS systems 
will form via a hypercritical common envelope phase. 

3. BIRTH SPINS OF BH IN BH-NS BINARIES 

No direct measurements of BH birth spins (or even 
BH spins) exist yet. A priori the BH birth spin cannot 
be constrained beyond the most fundamental level (i.e., 
to avoid naked singularities, we must have low angular 
momentum J: specifically, a = J/M^ < 1). For exam¬ 
ple, the collapse of a slightly-hypermassive NS rotating 
at breakup will, on dimensional grounds alone, produce 
a BH of spin of order half of maximal [a ^ 0.5] should it 
collapse; more realistic computations which allow for po¬ 
tential differential rotat ion support could push this value 
even higher (see, e.g.. iDuez et ^ 12004 iShaoiroi 120021: 
iHawke et a,1 .112110411 . 

Lacking unambiguous theoretical guidance for BH 
birth spins in merging BH-NS systems, and lacking direct 


observations of BH spin, we employ evidence for birth 
spin in similar systems to guide our choices for BH birth 
spin. On the one hand, on the basis of our understanding 
of NS birth spins and spinup in NS-NS binaries, we ex¬ 
pect low BH birth spins (a = 0). However, on the basis 
of suggestive observations of X-ray binaries, we suspect 
moderate BH birth spins could occur (a = 0.5) 

3.1. BH birth spin estimates from NS observations 

The same processes that produce BH birth spin (i.e. 
core collapse and accretion-induced spinup) also deter¬ 
mine the spins of the well-observed NS population. Ob¬ 
servations of the NS population can thus potentially pro¬ 
vide us with some constraints on BH birth spins. 

BH birth spin from core collapse: We expect we can 
estimate the birth spins a = J/M'^ of BHs, and particu¬ 
larly low-mass BHs, through the estimated birth angular 
momenta J of young NS. For BHs and NSs of compara¬ 
ble mass, the collapse process should be nearly identical; 
therefore, the collapse product should have very similar 
J and M. Assuming n o other process intervenes [for 
example, as discussed in lLindblom and Owed ll2002ll . r- 
mode damping is not expected to significantly change the 
NS angular momentum of young NS], NS will spin down 
electromagnetically, and thus NS spins at birth can be es¬ 
timated from the observed pulsar NS sample. In the past 
few years a number of studies of observed radi o pulsars 
have estimated the NS spin periods at bir th Isee lLorimeil 
iOOllK rame il l2003l: iMigliazzo et alJl2003. and references 
therein) in the range of 10 — 140 ms. These values are 
significantly slower than break-up spins (for typical NS 
equations of state) and correspond to a ~ 0.005 — 0.02 
(assuming a NS radius of 10km and rigid rotation). 

BH birth spin from accretion-induced collapse: Our 
simulations show that a significant fraction (i.e., 4199 out 
of 5435 binaries) of the BH in merging BH-NS systems 
were originally NS that experienced HCE and collapsed 
into BH. Therefore birth spins for this BH class are re¬ 
lated to the spin-up of NS during a CE phase. Our cur¬ 
rent understanding of this spin-up process is quite lim¬ 
ited. However, we can obtain guidance from the mildly 
recycled pulsars in known double neutron star binaries. 
These pulsars are believed to have been spun up dur¬ 
ing a CE phase (where collapse to a BH was avoided 
though) and it is evident that they are not spinning close 
to break-up: the fastest known pulsar in a double NS bi¬ 
nary is PSR J07 37-3039A spinning a t ~ 20 ms or else 
having a = 0.01 liBnrgav et a,lJl2f)0^ . Therefore post- 
CE NS spins appear to correspond to a no larger than 
a = 0.01. 

Based on the observational considerations discussed 
above, we would expect the birth BH spin to be negligi¬ 
ble (i.e., a < 0.03). However, a number of uncertainties 
in our modeling suggest that very different values of a 
could be equally plausible. 

Uncertainties in collapse model: For black holes that 
form from the supernova of a massive star, we assume 
that the birth spin of these BHs are comparable to the 
birth spin of comparable-mass NSs, and we then as¬ 
sume that those NS birth spins may be estimated from 
electromagnetic (pulsar) spindown from present-day pul¬ 
sars. However, most of these low-mass black holes form 
through fallback of post-supernova ejecta onto the hole; 
if the post-supernova fallback material carries significant 
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angular momentum and can spin up the nascent BH to 
high values of a, even if the protoneutron star is not 
spinning with a period shorter than ~ 10 ms. Further, if 
young protoneutron stars can spin down rapidly through 
other mechanisms, such as r-modes llTlndhlom and OwenI 
Enna, straightforward electromagnetic-based extrapola¬ 
tions may significantly underestimate the protoneutron 
star spin. If either of these mechanisms took place, then 
the birth BH spin could be significantly larger than the 
naive estimate we outlined above. 

Uncertainties in HCE model: The known pulsars in 
double NS binaries have avoided accretion-induced col¬ 
lapse into BH. Those that do not avoid this collapse must 
accrete a higher amount of mass and therefore it is still 
possible that they get spun-up to shorter spin periods 
and therefore lead to BH formed with higher a values 
than expected. 

Given these uncertainties, we cannot be sure whether 
the relevant BH birth spins are negligible (i.e., a = 0) 
or large, comparable to the break-up spin of a NS (i.e., 
a ~ 0.5). 

3.2. Black hole Spins in X-Ray Binaries 

Black hole X-ray binaries (XRBs) offer the possibility 
of directly measuring the spin of a compact object, since 
they involve a highly relativistic accretion flow in the 
strong field of a BH. Two techniques are prominent in 
the literature: interpreting quasiperiodic oscillations and 
fitting iron line profiles. 

Quasiperiodic oscillations are believed to be modes of 
the inner, highly relativistic regions of the BH accretion 
disk. As such, their frequency is expected to be related 
to the Kepler frequency of the inner edge of the disk, 
which in turn is intimately connected with the proper¬ 
ties of the black hole. Measurements of high-frequency 
QPOs that are inconsistent with the innermo st stable cir¬ 
cular o rbit of a nonspinning black hole, as bv IStrohmaveil 
1 I 2 OOIII and IHemillard et fflll2002ll . have been interpreted 
as evidence for BH spin. Moreover, evidence suggests 
the spin suggested by these measurements is not merely 
a relic of past accretion. Unfortunately, QPO frequen¬ 
cies cannot be unambiguously translated to black hole 
angular momenta; their interpretatio n depends strongly 
on disk mode mo delling (see, e.g. iTorok et al.l 120051: 
IRezolla et al.ll2?)?i^ . with results that vary from a « 0.1 
to a « 1 depending on the assumptions used. However, 
even if the QPO interpretation is correct, the associated 
BH spin estimates need not reflect BH birth spins; such 
spins are most probably reached via long-term disk ac¬ 
cretion, a process that is not relevant to BHs in BH-NS 
merging binaries. 

Iron line profiles potentially offer in principle a more 
direct probe into the inner disk rotation profile and, thus, 
the BH spin. The iron line can show strong Doppler shifts 
(due to rotation) and redshifts (due to strong gravity); 
in particular, very rapidly spinning BHs should show 
asymmetric red “tails”. While pioneered for use with 
AGN, thi s same technique has been a pplied to BH XRBs 
(see, e.sf.. IMiller et al.ll2002L I2fi04albfl . We treat the re¬ 
sults, several claims of extreme spins a > 0.8, with great 
caution however. On the one hand, these systems are 
strongly model dependent, and many physically relevant 
details (e. g., proper accounting of lig ht bending near the 
hole; see iBeckwith and Pond 120041 have yet to be in¬ 


cluded in the models. On the other hand, just like in the 
case of QPO interpretations, the estimated spins may be 
entirely due to the disk accretion in the XRB and there¬ 
fore unrelated to the birth BH spin. 


3.3. Birth Black Hole Spins Adopted 

To summarize, based on observations of isolated NS we 
believe the birth spins of NS - and thus of low-mass BHs 
- to be small. On the other hand, observations of X-ray 
binaries (QPOs and iron lines) suggests that higher-mass 
BHs (M > 5Mq) could be born with moderate spin. To 
allow for this substantial uncertainty, in what follows we 
will consider the implications of two models for the birth 
BH spin: a = 0 and a = 0.5. 


4. MASS ACCRETION AND SPINUP 


Figure ^ demonstrates that the black hole accretes a 
significant fraction of its mass during the HGE phase. 
Since that matter probably carries some amount of an¬ 
gular momentum, the black hole could conceivably spin 
up significantly during the HCE phase. Unfortunately, 
the details of HCE accretion ~ particularly at the fine 
level of detail needed to resolve the amount of angular 
momentum advected onto the compact object - are not 
understood. 

Lacking a quantitatively sound choice for the specific 
angular momentum, we argue instead that the specific 
angular momentum accreted should be at most the an¬ 
gular momentum of the marginally stable equatorial par¬ 
ticle orbit. This model, equivalent to assuming that the 
accretion proceeds through a thin disk, implies that the 
final black hole spin depends only on the total amount 
of matter accreted and the initial black hole spin J (or, 
as commonly deno ted, a = J/M^): as der ived initially 
bv iBarde^ iHQTflU (see also nr orn ilMi , for a black 
hole which is initially nonspinning with mass Mi, the 
BH spin parameter a is given by the following expression 
a = aB{M,M,y. 


aB{M,Mi) 



4- 




• ( 3 ) 


We consider this estimate to represent an upper limit to 
the BH a, since we do not expect that accretion in a HGE 
eyent will occur through a thin disk. Instead, we expect 
that on almost all scales the accretion flow will be nearly 
radial. 


4.1. Accretion onto nonspinning holes 

We use our results on compact object masses at the 
onset and the end of HCE to estimate a under the thin 
disk assumption. We consider two cases: (i) all mass 
accreted during the CE phases contributes to spinup and 
(ii) if accretion-induced collapse (AIC) occurs, only the 
fraction of the mass accreted after AIC contributes to 
spinup. 

All mass contributes to spinup: As demonstrated in the 
top panel in Fig. [3 because of the substantial amount of 
mass accreted during HCE phases, nearly all BHs could 
conceivably spin up to very large a; a histogram of the 
values of a we observe would be highly concentrated near 
a ~ I. In particular, many low-mass BHs (i.e, those 
which formed via AIC during an HCE phase) have ac¬ 
creted of order half their final mass, and can potentially 
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Fig. 2. — For systems that undergo HCE evolution (i.e., 
the “standard” or “reversed” channels presented in Table [Q this 
plot shows the BH spin expected assuming the HCE accretion 
occurs through a thin equatorial corotating disk: points are 
(Ma final? a(AM)) where a is computed using the Bardeen formula 
[Eq. m assuming an initially nonrotating hole. The two panels 
show the results assuming (i) all the mass that is accreted con¬ 
tributes to spinup and (ii) only the mass accreted after AIC of 
the NS into a BH. The solid curves (shown only for M > 4) show 
the empirical upper bound presented in Eq. 11) for nonspinning 
newborn BH. 


spin up to of order a ~ 1. Further, every system that 
undergoes HCE evolution accretes enough mass to spin 
up significantly (a > 0.2). In this case, the BH-NS bina¬ 
ries with the largest spins (i.e., which have a ci 1) have 
mass ratios between 2 : 1 to 3.5 : 1. 

However, for NS which undergo AIC, this approach 
approach will strongly overestimate the maximum at¬ 
tainable spin, since it assumes disk accretion at the 
schwarzchild radius, even before the BH forms. 

Only post-AIC mass contributes to spinup: Since cur¬ 
rent astrophysical evidence suggests that NS do not get 
strongly spun-up during CE events, we also consider the 
possibility that only the mass accreted after the BH is 
formed through AIC contributes to its spin-up. (Those 
BH which form through direct collapse are treated as in 
the previous case.) In this case, summarized in the bot¬ 
tom panel in Fig.|21 we find a much broader distribution 
of possible a values (i.e., a histogram of possible a out¬ 
comes would be more nearly flat, from a = 0 to a = 1). 
In this case, the systems with the largest spins have mass 
ratio near 3.5 : 1. 

Both models have the same behavior for systems with 
large final BH masses (Mbh,f > dM©). Good empirical 
lower and upper bounds for a for Mbh,f > 4 Mq fo llows 
directly from the mass-accretion bounds Eqs. (I1I2II : 

a<aB^Mhhj,—^^ — 8 Mq/5^ (4) 

a>aBiMthj,0.95Mbhj-0.57MQ) (5) 

where and Mbh,i are the BH final and birth masses, 
respectively. 


Fig. 3. — As Fig. 121 but on a different vertical scale and assuming 
an initially rotating hole: a = 0.5. The solid curves (shown only 
for M > 4) show the empirical upper bound presented in Eq. m- 


4.2. Accretion onto spinning holes 

Black holes which are initially spinning present a more 
complex accretion challenge in principle: the BH spin 
need not be aligned with the disk angular momentum 
axis, so the accreted material could just as well spin down 
as spin up the hole. However, following the same process 
we used to circumvent the the (substantial) uncertainties 
we addressed earlier, we limit attention to bounding the 
BH spin. The most conservative bound is obtained by 
assuming a corotating equatorial disk. 

Bardeen’s formula [Eq. Q] applies equally well to 
BHs which are initially spinning, if Mi is chosen so the 
true initial black hole mass M^h^i and spin Oi satisfy 
aB{Miih,i) = cLi. To be concrete, if the birth BH spin 
is Gi = 0.5, then choosing 

M, = 0.84 Mbh,^ (6) 

yields aB{Mbh,i,Mi) = 0.5. When we repeat the analysis 
of the previous section, we find results summarized in 
Figure El Since Bardeen’s spinup relation is monotonic, 
the empirical upper limit presented in Eq. 0) [for Gi = 0] 
translates directly into a corresponding spin limit when 
Gi = 0.5 


a < as j, 0.84 x 

G>aB (Mbhj, 0.84 X [0.95Mthj - O.STMg]) (8) 

Briefly, the results are qualitatively similar to the non¬ 
spinning case, but compressed in scale such that the min¬ 
imum possible spin is a = 0.5 instead of a = 0. Since the 
construction of spin as a function of BH mass is simi¬ 
lar, the spin upper bound also follows Eq. 0 , with the 
Bardeen formula calculated choosing Mi so the initial 
black hole has spin a = 0.5. 


'^Mbhj 


- 8Mq/5 


■ (7) 


4.3. General Conclusions 
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Our simulations have shown that almost all merging 
BH-NS binaries form through a hypercritical common 
envelope (HCE) phase, where they accrete a substan¬ 
tial fraction of their final mass. Lacking a quantitative 
model for the specific angular momentum accreted dur¬ 
ing HCE, we use the most pessimistic model: we assume 
the BH spins up through a thin equitorial disk, accord¬ 
ing to the Bardeen formula [Eq. ||21)]. Even setting aside 
those systems with even greater physical uncertainties 
- those systems with final BH mass M < 4:Mq, which 
form through accretion-induced collapse of the NS and 
whose final spin depends additionally on the details of 
the collapse and angular momentum transport onto the 
NS during that collapse - and even ignoring the possi¬ 
bility of BH birth spin, we find that large spins (a ~ 1) 
could easily be attained in principle from the HCE phase 
alone: with such a substantial fraction of its final mass 
transferred, the BH could easily spin up. 

However, for BHs with M > dM©, significant spinup 
after birth becomes increasingly difficult, as larger black 
holes require ever larger (and more unlikely) mass trans¬ 
fers in order to spin up to large a. For example, assum¬ 
ing nonrotating initial BHs, the only BHs we found with 
spin a > 0.8 are those with final masses between 3Mq 
and IOMq; similarly, for BHs with a = 0.5 at birth, the 
only BHs we found with spin a > 0.95 are those with 
final masses in this same range. 

5. BH SPIN AND DETECTION OF BH-NS INSPIRAL 

In BH-NS binaries, the spin of the BH can, 
if large, have a significant effect on the late 
stages of inspiral and in particular the e r nitted 
gravitational waveforms fsee lAn ostolato s"^ 
ApostolatojllOObt iGrandclement et al.l 12004 iPan et al.1 
2004 iBuonanno et al.ll20(M and references therein). As 
a result, in BH-NS binaries measurements of gravita¬ 
tional waves offer the potential to extract information 
about the spin. Furthermore, because of precession, such 
binaries will be much more reliably detectable, suffering 
less frequently from the effects of poor source orientation. 

But these advantages come at a price. The more 
complex waveforms of spinning binaries are vastly more 
complicated to model. At present only ad -hoc mod¬ 
els a re considered computationally feasible l|Pa,n et af.l 
l2f)f)l , and even these methods involve substantial com¬ 
putational challenges, particularly at more extreme mass 
ratios. 

Our calculations provide extremely conservative upper 
bounds on the BH spin in BH-NS binaries. These bounds 
tell us that employing spin need not pose such a dramatic 
dilemma, on astrophysical grounds: if we employ spin, 
we could require fewer computational resources than we 
might expect; and if we neglect it, we will not lose a 
dramatic fraction of potentially observable events. 

5.1. Mass ratio spin constraints and detection 

When constructing template banks for BH-NS inspiral 
a priori, parameter ranges are typically chosen to cover 
a very conservative parameter range; for example, in one 
of the few papers that p resented results for variable mass 
ratio, IPan et ;il] ll2004l) built a template bank based on 
all binaries with masses in the range Mhh S [7, 12]Mq 
and Mns & [1, 3]Mq and BH spins in the range a G [0,1]. 

However, most of the templates in the template bank 


come from extreme-mass-ratio systems, involving a high- 
mass BH which we have just demonstrated is difficult to 
spin up. Therefore, a different template bank, which ne¬ 
glected astrophysically irrelevant systems, could perform 
the same search much more rapidly. Detailed computa¬ 
tions of the computational savings we could obtain by 
moderating our template bank are far beyond the scope 
of this paper, and will depend sensitively on two poorly 
constrained parameters, the BH birth spin and the max¬ 
imum NS mass (here, assumed M = 2Mq.) However, 
in an optimistic case we expect that a more judicious 
choice of templates could reduce the template bank size 
by a significant factor, cf. Sec. VI C of Pan et al. 

5.2. Limits on the decrease in event rate due to 
omitting BH spin 

For the purposes of detection, spinning templates are 
required only when nonspinning templates fail to have 
sufficient overlap with the physical signal to guarantee 
most inspiral events are detected. Loosely speaking, a 
nonspinning template fails to mimic the effect of preces¬ 
sion on the waveform. Because of spin-orbit coupling, 
angular momentum will be exchanged, causing the or¬ 
bital plane to process around the total angular momen¬ 
tum. Therefore, spinning templates are needed to de¬ 
tect those systems where the effects of precession are 
strongest - namely, those with (i) strong BH spins awl 
that are (ii) strongly misaligned wi th the orbital angu¬ 
lar m omentum (see, e.g.. Fig. 2 in IGrandclement et alJ 

l2^. _ 

IGrandclement et ~ ] il2nnl have quantitatively exam¬ 
ined the expected overlap between nonspinning tem¬ 
plates and the expected signal from binaries in which the 
BH spin magnitude and orientation are independently 
varied (see their Figure 2), assuming BH-NS binaries 
with Mbh = IOMq and M„s = 1.4Mq. They then fur¬ 
ther convolved this function with an expected spin orien¬ 
tation distribution to derive an estimate for the fraction 
of events that could be seen that will be seen, as a func¬ 
tion of spin magnitude of the BH (see their Figure 3)."^ 
They find that even if all BH were spinning maximally, 
only about 30% of potentially-detectable BH-NS mergers 
will be missed. 

The Grandclement et al. results, however, are not di¬ 
rectly applicable in our circumstances. Spin-dependent 
modulations depend significantly on mass ratio, a param¬ 
eter that is varied in our BH-NS binaries but fixed for 
those of Grandclement et al. Nonetheless, if we employ 
their relation between (i) the probability of detecting a 
detectable source when ignoring spin and (ii) the BH spin 
magnitude of that source, then we find when we convolve 
with our Monte Carlo sample of BH spins - which, re¬ 
call, was designed to produce the largest plausible spin 
values, given known mass transfer - that we would lose 
at most 30% of events in extreme cases (i.e., high BH 
birth spin or all HCE mass accreted onto NS) and slightly 
less (25%) in more conservative cases (i.e., mass before 
AIC does not contribute to spin). Briefly, searches which 
use nonspinning template banks will likely find most de¬ 
tectable inspiral signals. Practically speaking, given the 
considerable uncertainty in the underlying event rate it- 

^ Specifically, they plot the average of the cube of the fitting 
factor versus the BH spin a. 
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self, upper limits produced using precise templates will 
have only marginally more astrophysical impact than up¬ 
per limits produced using nonspinning template banks. 
Of course, for parameter estimation, correct templates 
are essential in order to extract the BH spin, and thus to 
constrain its accretion history. 

6. SUMMARY AND FUTURE DIRECTIONS 

In this paper we perform population synthesis calcu¬ 
lations to track the history of systems which evolve into 
merging BH-NS binaries, demonstrating that these ob¬ 
jects form through only a few channels, almost exclu¬ 
sively involving a hypercritical common envelope (HCE) 
phase and often involving accretion induced collapse of a 
NS into a BH. We show that a significant amount of mass 
is typically accreted during this HCE phase. Using an 
upper limit based on thin-disk aecretion, we find the BHs 
accrete enough matter to potentially spin up to large a 
values, independent of the BH birth spin. Finally, we 
note the observed sample of neutron stars suggests the 
birth spins of the black holes in these binaries should be 
small; therefore, the BH spins in BH-NS binaries should 
arise almost entirely from the matter they accrete in a 
HCE phase. Unfortunately, the HCE phase is poorly un¬ 
derstood. While our result is a first step of a program to 
better understand the expected BH spins for those sys¬ 
tems LIGO and other gravitational-wave detectors could 
observe, more simulations of these processes are needed 
to obtain firmer conclusions. 

Further, based on our rough understanding of the 
range of expected spin magnitude, and using the ex- 
pected effect of sp i n on GW detection drawn from 
iCrandclement et alJ ll2fln4fl . we demonstrated that ne¬ 
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glecting spin in GW searches should only moderately 
reduce the detection rate. However, our analy sis was 
based on the work of iGrandclement et alJ ll2004ri . which 
like many papers in the gravitational-wave literature as¬ 
sumed all BHs in BH-NS binaries had M = IOMq. In 
fact, based on our calculations we expect (i) the mass ra¬ 
tio can vary, and should be biased towards larger values 
(i.e., closer to 3 : 1 than 10 : 1.4); and (ii) the sam¬ 
ple of BH-NS binaries should show strong correlations 
between mass ratio and spin magnitude, because higher- 
mass BHs are harder to spin up. We hope to undertake a 
more thorough Monte Carlo study of the effect of various 
expected distributions for BH spin magnitude, spin tilt, 
and BH-NS mass ratio in a future paper. 

Finally we note that in our simulations we consider 
only BH-NS in the Galactic field; we do not account for 
any stellar interactions relevant to centers of globular 
clusters and similar dense environments. If BH-NS bina¬ 
ries form in significant numbers in cluster centers, then 
their spin properties could be entirely unrelated to the 
analysis presented in this study. 
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